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ABSTRACT 


The  MGR[v]  algorithm  of  Ries,  Trottenberg  and  Winter,  Algorithm  2.1 
of  Braess  and  Algorithm  4.1  of  Verflirth  are  all  algorithms  for  the  numerical 
solution  of  the  discrete  Poisson  equation  based  on  red-black  Gauss-Seidel 
smoothing  iterations.  In  this  work  we  consider, the  extension  of  the  MGR[0] 
method  to  the  general  diffusion  equation  -V  •  pVu  =  f  v.  In  particular,  for 
the  three  grid  scheme  we  extend  ,an  interesting  and  important  result  of  Ries, 
Trottenberg  and  Winter  whose  results  are  based  on  Fourier  analysis  and  hence 

n  ■  ■ '  »  :  '  't 

intrinsically  limited  to  the  case  where  ft  is  a  rectangle.  Let  ft  be  a 

t  •  r  -  .  "  ■ 

general  polygonal  domain  whose  sides  have  slope  f 1 ,  0  and  <p  :  Let  eu 

be  the  error  before  a  single  multigrid  cycle  and  l.et-e  be  the  error  after 

this  cycle.  Then  ||  e1 1| .  £  \  (1+Kh)  ||  e°||  where  ||  || .  denotes  the 

Lh  Lh  Lh 

energy  or  operator  norm.  When  p(x,y)  =  constant,  then  K  =  0  . 


.• ..  :  t  I 


MOV  1  3  1985  ■  J 


t 
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1 .  Introduction 

The  MGR[v]  multigrid  algorithms  of  Ries,  Trottenberg  and  Winter  [4], 
the  Algorithm  2.1  of  Braess  [l],  [2]  and  Algorithm  4.1  of  Verflirth  [5] 
are  all  algorithms  for  the  numerical  solution  of  the  discrete  Poisson 
equation  (the  usual  5-point  difference  equations  with  Ax  =  Ay  =  h)  based 
on  red-black  Gauss-Seidel  smoothing  iterations.  The  analysis  of  [4]  is 
based  on  Fourier  Analysis  and  is  restricted  to  the  case  where  the  basic 
domain  ft  is  a  square.  The  analysis  of  [1],  [2]  and  [5]  is  for  a  bounded 
polygonal  domain  ft  whose  sides  have  slope  +  1,0  and  °°  and  is  based  on 
certain  energy  estimates  and  a  particular  interpretation  of  the  matrix  equa¬ 
tions.  While  this  is  not  explicitly  stated,  this  interpretation  can  be 
viewed  as  a  particular  choice  of  I ^ ,  1^ ,  1^,  I^h  etc,  the  operators  which 
carry  on  the  communication  between  the  grids. 

Recently,  Kamowitz  and  Parter  [3]  considered  a  generalization  of  the 
algorithms  of  Ries,  Trottenberg  and  Winter  and  Braess.  They  consider  the 
general  diffusion  equation 

V  0  p(x,y)Vu  =  f  in  ft  , 

^  1 )  u  =  0  on  9ft 

p(x,y)  >  P0  >  0  , 

in  general  domains  ft  .  Using  a  different  choice  of  l|j,  1^  than  Braess, 
i.e.  imagining  a  different  interpolation  structure  in  the  space  S^,  they 
employ  other  "Energy  Estimates"  to  obtain  the  basic  estimate  -  for  a  two 
grid  scheme:  let  denote  the  error  before  a  single  multigrid  cycle  and 
let  denote  the  error  after  that  complete  multigrid  cycle,  then 


2 


0.2)  II  £\  <  i-tl+Kh)  ||  e°||L 

h  h 

where  the  constant  K  depends  only  on  Pg  and  HVpU^,  the  ®  norm 

of  the  gradient  of  p(x,y)  and  ||  || ,  denotes  the  operator  or  energy 

Lh 

norm.  However,  it  is  important  to  remark  that  despite  the  different 
interpretation  of  the  problem,  in  the  case  of  constant  diffusion  coefficient 
p(x,y)  =  1  we  are  dealing  with  exactly  the  same  problem  and  the  same  itera¬ 
tive  method.  The  estimate  (1.2)  is  thus  a  generalization  of  the  estimates 

(1.3)  p(MG)<£,  p(MGR[0])=y 
of  [l]  and  [4]. 

Another  remarkable  estimate  of  Ries,  Trottenberg  and  Winter  [4]  is  the 
fact  that,  in  the  case  of  Poisson  equation  in  the  square,  if  a  third  grid 
is  introduced  and  one  uses  the  MGR[0]  method  one  obtains 

(1.4)  p(MGR[0] ,  3  grid)  =  \  . 

In  this  report  we  obtain  this  estimate  in  the  form  (1.2)  for  the  general  dlf 
fusion  equation  (1.1)  in  bounded  polygonal  domains  ft  whose  sides  have  slope 
+1,0  or  <*.  We  also  require  that  the  corners  of  ft  belong  to  the  coarsest 
mesh.  The  constant  K  is  a  constant  depending  only  on  Pg,  and  the  °°  norm 
of  the  first  and  second  derivatives  of  p(x,y).  Moreoever,  if  p(x,y)  =  const, 
then  K  =  0  .  In  general,  throughout  this  paper  K  will  denote  such  a  constant 
In  section  2  we  formulate  the  problem  and  the  basic  three-grid  multigrid 

A 

iteration.  In  particular  we  introduce  the  coarse  grid  operators  L^,  L^,  l^, 

A 

I^-  In  section  3  we  develope  more  notation  and  recall  some  basic  estimates 
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of  [3],  In  this  section  the  reader  is  introduced  to  a  number  of  additional 

difference  operators  ,  L^,  L2J^,  L^,  Qx»  Mx,  Lx  .  This  plethora 

of  operators  gets  a  bit  confusing.  However  if  one  first  concentrates  on  the 

case  p (x,y)  =  1  (i.e.,  the  Poisson  equation)  the  situation  simplifies.  In 

this  case  ,  l_2h  =  and  (we  always  have)  ^ , 

-  i  (i)  wn 

*"2h  =  2zh  +  2"L2h  '  Moreoever,  in  this  case 


L  =  Q 
x  ^x 


L  =L(1) 
L2h  lH 


q 


2h 


C^2h  is  the  coarsest  grid]  and 


Vfl2h 

is  the  intermediate  grid].  Another  observation  which  should  be  useful 
is  the  fact  that,  in  this  case  is  the  same  difference  operator  as 

L2h  excePt  f°r  P01'nts  in  which  are  next  to  the  boundary.  Moreover, 
these  exceptional  points  are  in  nH/^2h  not  in  ^h’  This  perturbation  of 
causes  a  technical  difficulty  in  the  proof  of  lemma  5.2  even  in  this 
simplest  case.  In  all  cases  the  introduction  of  the  variable  diffusion 
coefficient  p(x,y)  introduces  perturbation  of  the  basic  operators.  However, 
the  essence  of  the  proof  of  the  main  result  [Theorem  5.1  or  the  estimate 
(1.2)]  is  contained  in  the  constant  coefficient  case.  The  analysis  of  the 
algoritfm  is  given  in  two  parts,  sections  4  and  5. 
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2.  The  Problem 

Given  a  (small)  value  h  >  0  let  {(x.  ,y.)  =(kh,jh);  k,j  =  0,+l ,+2, . . . } 

K  J 

be  the  associated  mesh  points  in  the  x  -y  plane.  Let 

(2.1)  Rq  :  =  {(xk^j);  k+j  =  l(mod  2)} 

(2.2)  Rg  :=  {(xk.yj);  k  =  j  =  0  (mod  2)} 

(2.3)  Rg  :=  {(xk,yj.);  k  5  j  5  ]  (mod  2)}  . 

Let  ft  be  a  bounded  polygonal  domain  in  the  plane  whose  sides  have 
slope  +1,0,  or  °°,  and  every  corner  point  (x,y)  of  9ft  belongs  to  R_. 

D 

Defi ne 

(2.4a)  fth  ■  (RQURBURe)  O  ft 

(2.4b)  9fth  =  (R0URgURG)  n  3ft 

(2.5a)  ftH  «  (RbuRg)  n  « 

(2.5b)  9ftH  =  (Rg  uRq)  n  3ft 

(2.6a)  ft^^  =  Rg  Hft 

(2.6b)  9ft.,,  =  R  naft  . 

tn  B 

For  any  function  F(x,y)  defined  on  ft  we  write: 
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(2.7a) 

Fk,j  "  F(xk,yj)  • 

(2.7b) 

Fk*.J  *  F((kt!i)h'yj>  • 

(2.7c) 

Fk,M  *  F<V'J^>h>  • 

The  algebraic  problem  to  be  solved  is:  Find  a  mesh  function  U  =  (U^ 

defined  on  ft.  u  3ft.  which  satisfies 
h  h 


(2.8a)  [L„U]kJ  =  Fkj  ,  Uk>,.)  £  Oh 

(2.8b)  Ukj  =  0  ,  (xk,yj)  €  3fth 


where 


(2.8c) 


[LhU]kj 


=  h2  {pk-%,j^Uk,j  _Uk-1,jl  '  pk+%,  j^Uk+l ,  j 
^2  ^Pk,  j-^Uk,  j  "Uk,j-1-'  '  Pk, j+J}^Uk, j+1 


■Uk,j 


Uk,j 


]}  + 

]}  . 


We  turn  to  solution  of  these  linear  algebraic  equations  by  a  three- 
grid  method. 

Let  S^,  S^,  S^h  be  the  linear  spaces  of  mesh  functions  defined  on 
fth  u  3ft^,  ft^  u  3  2^  and  ft^  u  sft^^  respectively  which  vanish  on  the 
respective  boundaries  3ft^,  Sft^,  3^2^  •  We  set  up  communication  between 
these  spaces.  Specifically  we  define  the  linear  interpolation  and  pro¬ 
jection  operators  ijj,  I^,  I ^ ,  1^  as  follows.  The  interpolation 
operator  1^  (see  the  definition  of  of  [3])  is  given  by 


(2.9b) 

Wtj  -  V  <f(xk.yj)^Hu 

3Q^I 

and,  if 

(x^,y.)  £  »  then 

(2.9c) 

[IHU]kj  ck.  {pk-is,jUk-l,j  +Pk+is,jUk+l,j 

+ 

Pk,j-hUk,j-l  +Pk,j+i5Uk,j+l} 

where 

(2.9d) 

Ckj  "  {pk+ij,j  +pk-ij,j  +Pk,j-J5  +Pk,j+J5} 

• 

Of  course,  if  (x^y^. )  e  3^/3^  then 

(2.9e) 

C'Kj  - 0  • 

|_| 

The  projection  operator  1^  is  defined  by 

(2.10) 

IH  =  l(ih)T 

ih  2uH;  • 

Remark:  The  factor  Jr  in  (2.10)  is  included  merely 

consistent  with  the  MGR[v]  methods  of  [4]. 

to  keep  the  method 

The 

f_J 

interpolation  operator  I?.  is  defined  in  a 

similar  manner  by 

7 

with 

(2.11b) 

^2hU^kj  =  Ukj  1  if  (xk,yj)  e  Q2hU3Q2h  * 

and,  if 

( ,y j )  €  ^/^2h  *  *  *  then 

(2.11c) 

H  1 

*-I2hU-^kj  ck.  tpk+is,j+i5Uk+l  ,j+l  +  pk+Js,j-JsUk+l  ,j-l 

+  Pk-i5,j+i5Uk-l  ,j+l  +  Pk-4,  j-JsUk-l  ,j-l  ^ 

where 

:  (2. lid) 

^kj  ”  ^pk+ig, j+^s  +Pk+^,j-h  +Pk-J5,j+i5  +Pk-is 

and,  if 

(x^.yj)  e  >  then 

(2. lie) 

[‘>kJ  ■  0  • 

i  The  projection  operator  I2h  is  given  by  1 

:  (2.12) 

I2h  =  1(IH  )T 
iH  2  u2h;  . 

i  Finally  we  define  the  "coarse  grid"  operators  L^,  L,,^  .  These  are 

(2.13a) 

LH:  SH  SH 

,  where. 

if  (xk,yj)  e 

■  (2.13b) 

^LHU^kj  “  2h2  {ck,juk,j  Pk+H,j+J5Uk+l  ,j+l 

i 

'  Pk+*5,  j-^k+l  ,j-l  "  Pk-15,j+i5Uk-l ,  j+1  ~  pk-*5,j-*sUk-l  ,j-l } 

► 
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and 


(2.14a) 


L2h:  S2h  "  S2h 


where,  if  (x^.y.)  e  then 

(2.14b)  [L2hU]kj  =  -V  tpk_,  f[U 


4h 


I  lpk-l,jLuk,j  -Uk-2,j]  '  pk+l,j[lW,j  'Uk,j]) 


*  ^2  {pk,j-ltUk,j  -Uk,j-2^  'pk,j+l[Uk,j+2  'uk,J]>  ' 

We  are  now  ready  to  describe  the  three  grid  methods.  Let  Bh  be  a 
non-singular  linear  operator  defined  on 


(2.15) 


Bh:  Sh  ^  Sh  - 


Let  the  smoothing  operator  Gh  be  defined  by 


(2.16a) 


G„  -  !„  -  B''lh 


and  assume  that 


(2.16b) 


<LhGhu>Ghu  } 
< Lhu,u  > 


<  1  , 


Vu  €  Sh,  u  M  . 


Al qori thm 


Step  1  : 


Given 


form 


~u  =  Ghu 


0 


F  . 


(2.17) 
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3.  Some  Notation  and  Facts 

Let  u,v  e  or  or  . 


Then 


(3.1)  l  “k./k.j 

where  the  sum  is  taken  over  all  indices  (k,j)  so  that  (x^.y^.)  e  ,  or 
or  ^2h  resPectl‘ve1y-  Whenever  it  seems  that  further  clarity  is  required 
we  will  indicate  the  space  by  writing 


<  u,v  >  ,  a  =  h  or  H  or  2h  . 

a 

Since  L^,  and  L£^  are  positive  definite  operators  we  have  the 
inner  products 


(3.2) 

[u.v], 

< 

<Lau'v>a  ' 

a  =  h  or  H  or 

Let 

(3.3a) 

V 

=  Null  space 

T"|_  C  S 
h  h  ^h 

(3.3b) 

Kh 

:  =  Range  l|j 

csh 

(3.3c) 

Nh: 

=  Nullspace 

ll\  C  SH 

(3.3d) 

kh 

:  *  Range  l”h 

CSH 

Lemma  3. 

]_:  We  have 

(3.4a) 

Sh 

■  Nh  ©  Rn  ■ 

s„  •  nh©rh. 

2h  . 


In  fact,  N 

h  and 

]R^  are  L 

^  orthogonal 

;  Nh  and 

orthogonal . 

That 

is,  if 

n  e 

N  ,  w  e  JR  , 

a  a 

a  =  h  or 

(3.4b) 

Cn 

,u]a 

’  aan-“>a  = 

0  . 

A  function 

<1  £ 

is  in 

Rh 

if  and  only 

if 

(3.5a) 

kj 

0,  (xk,yj) 

A  function 

V‘  Sh 

is  in 

Nh 

if  and  only 

if 

(3.5b) 

VkJ  = 

0,  (xk,yj). 

£  . 

A  function 

U  £  SH 

is  in 

if  and  only 

if 

(3.6a) 

CLHU]k 

0.  <vV  1 

£  ftH/fi2h  . 

A  function 

V  £  SH 

is  in 

nh 

if  and  only 

if 

(3.6b) 

Vkj  = 

0 , 

(xk,y^)  e  n 

2h  * 

are  L„ 
then 


Proof:  The  assertions  (3.5a)  and  (3.6a)  follow  from  the  definition  of 
l|j,  I2h  etc.  given  by  (2.9)— (2.1 2) .  The  assertions  (3.4a),  (3.4b),  (3.5b), 
(3.6b)  now  follow  immediately.  See  [3].  I 


Let 


(3.7a) 

A 

lh: 

=  l\  I 

Vh1 

(3.7b) 

A 

L2h: 

= 

Using  the  basic  relations  (2.10),  (2.12)  we  see  that 
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(3.7c) 


IHV II  L.  =  <LhIHv'IHv>h  ’ 
h 


(3.7d) 


1^*  LI  II  ^  =<LI^UI^U)  =  —  <  I  0  0) 
i2hU|lLH  LH12hu,i2hu  H  2  L2hu,u  2h  ’ 


The  formulae  (2.9),  (2.10),  (2.11),  and  (2.12)  together  with  (3.5a)  and 
(3.6a)  imply 


(3.8a) 

(3.8b) 


P  1  ,  A. 

LHU  "  2  Lh!H 


ft. 


4hv  =  2  LHT2hv 


ft 


2h 


A  /\ 


The  analysis  of  [3]  is  based  on  the  following  facts  about  L^,  l_2h  . 
Lemma  3,2:  There  are  operators  ,  L^,  such  that; 


(3.9a) 


lH  2  lH  2  lH  ’ 


(3.9b) 


L2h  2  L2h  2  L2h  ’ 


The  operator  is  based  on  the  five  points  (x^y^.),  (xk+^ ,yj+1 ) , 

(xk_i  ,yj+i ) ,  (xk_i  )  >  (xk+1 » y j _ i )  •  These  are  the  same  points  on  which 
Lh  is  based.  The  operator  is  based  on  the  five  points  (x^y^), 

(xk+2’yj)’  ^xk-2,yj^’  ^xk’yj+2^’  ^xk,yj-2^'  If  k  =  J  5  0  (mod  2),  these  are 
the  same  points  on  which  l2h  is  based.  Similarly,  if  k  =  j  e  0  (mod  2), 
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(j) 


/ 


\ 


tk) 


tk) 


The  five  point  star  for  L^.L^ 


0) 


(j) 


(k) 


■  denotes  a  point  in  Rg 
□  denotes  a  point  in  Rg 


The  five  point  star  for 


i 

Figure  1 


I 
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Is  based  on  these  same  points.  The  operators  ,  L ^  are 

"almost"  the  operators  L^,  l^.  To  be  precise,  we  have:  let 


(3.10a) 


Ms.j-ij 


Pv_U  iP 


k-»i,JKk-l,J-is  +  Hc.j-^k-^j-l 


'k-1  ,j 


'  k ,  j  - 1 


(3.10b) 


k+%,  j-*s 


Vj-^k-HsJ-l  +  PkjfrJpk+l,J->s 


'k, j-1 


k+1  ,j 


dkj  [ak-h,j-h  +  a  k  +^s ,  j  +  bk+Js,j-35  +  bk-*5,j+^ 


If  (k+j)  e  0  (mod  2),  then 

(3.11)  tLH1)u^kj  3  72  {"ak+is,j-H5Uk+l,j+l  "  ak-35,j-S5Uk-l  ,j-l 


-  bk+Js,j-35Uk+l,j-l  ■  bM5,j+J5Uk-l,j  +  l  +  dkjUkj} 


An  easy  computation  shows  that 

‘  -  Kh2 

'Vh-Vh1^'  • 

Hence,  for  every  U  e  S^, 

(3.12a)  |<LhU,U>  -  <lJ1}U,U>  |  <  Kh2  <LHU,U>  , 

( <  LhU,U>  -  <L^Vu>  I  <  Kh2  <lJ])U,U>  . 


(3.12b) 


A  basic  estimate  is:  for  every  U  £  S^, 


(3.13) 


0<(lJ1Vu>  1  2(1  +Kh)  (L^U.U) 


Hence,  if  we  write 


(3.14a) 


[  .  1l 

H  2  H  2  H  ’ 


then 


(3.14b)  -Kh<!_HU,U>  <  <L^U,U>  <  2(1+Kh)  <  LhU,U> 


Similarly,  let 


(3.15a)  A 


k+1  ,j 


Pk+*s,j+%pk+3/2,j+>s  +  pk+*s,j->sPk-3/2,j-*s 


'k+1 ,  j+1 


'k+1  ,j-l 


(3.15b)  B 


k,  j+1 


pk+^,j+^pk+^,j  +  3/2  +  pk->s,j+*spk->s,j+3/2 


'k+1 ,  j+1 


'k-1 ,j+l 


(3.15c)  Dkj  =  CAk+1>j  +  Ak_ljj  +  BkJ+1  +  Bk>j._1J 


If,  k  =  j  =  0  (mod  2) , 


(3'16)  L2h}  =  2^2  {"Ak+1 ,Uk+2,j  ‘  Ak-1 , jUk-2, j  ‘  Bk,j+lUl 


'  Bk,j-lUk,j-2  +  Dk,jUk„ 


>  j+2 
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An  easy  calculation  shows  that 

<3-,7a>  |2Aw,rW-Kh2> 

(3J7b)  |2Bkijtl-PkiJ+,liKh2. 

Hence,  for  all  U  e  S2^ 

(3.17c)  l<L2h)u’U>2h  “  <L2hU*U>2hl  -  Kh2  <L2hUsU>2h  • 

The  anolog  of  the  basic  estimate  (3.13)  holds.  That  is 
(3.18)  0  <  <l£JVu>  <  2(1+Kh)  <L^U,U>  . 

Hence,  if  we  write 

(3.19a)  4h  =  1  L2h  t  i  L<2> 

then 

(3.19b)  (-Kh)  <L2hU,U>  <  <L^U,U>  <  2(1+Kh)<  L^U.U  >  . 

Of  course,  if  p(x,y)  =  1  ,  then 

(3  20)  L  =  ^  L  =  ^ 

Proof:  The  construction  of  and  the  basic  estimate  (3.13)  is  found 

in  [3].  The  construction  of  and  the  estimate  (3.18)  then  follows 

from  the  same  arguments.  The  estimates  (3.11),  (3.17)  are  direct  computa- 
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Our  next  result  looks  at  the  operator 


(1) 


Lemma  3.3:  The  operator  is  of  the  form 


(3.21)  [^]kj  =  g  {-\+l,jUk+2,j  ‘  \-l,jUk-2,j  '  ^k,j+lUk, j+2 


-  B,  •  iU.  .  ,  -  D.  .11.  . 

k »J-1  k,j-2  k,j  k,jj 


The  coefficients,  A,  B,  D  are  given  by 


(3.22a)  A 


k+1  ,j 


1  2 

_  Pk+*s,jpk+3/2,j  +  2  ^Pk+H,j^  9 k+1  ,j 


Ck+1 ,  j 


(3.22b)  A 


k-1  ,j 


(3.22c) 


'k-1 ,  j 


1  2 

Pt,  4-U.Pl,  4*3,.  +2(Pk. -j+J  0k..i- 


(3.22d) 


5  =  Hc..HV'k,J  +  %  +  2  ^k.j+V  k,  j+1 

k>J+1  ck, j+1 


B  =  Pk,j-)sPk,j-3/2  *  2(pk,j-^  6k,j-l 
k,j-l  c, 


'k,  j-1 


(3.22e)  =  Ak+1>j  +  ^  ^  +  §kJ_]  +  §kJ+] 


where 

(3.23) 

Proof: 

Remark: 

Lemma  3. 

(Vi,yj 

Proof: 
that  all 

It  i 
based  on 
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e  =  \ 

U,a 


0,  (xM,y0)  i  3Sih 


These  coefficients  were  computed  in  [3] 


If 


k±l ,  j  *  0  ’ 

then  Uk±2,j  =  0  ’ 

k,j±1  ^  0  * 

then  Uk  .12  -  0  . 

4:  Let  (x^.y^)  e  ^  •  Then  all  4  of  its  h  grid  neighbors 
) .  (xk,yj±1 )  e  .  Hence 


ek±l,j  =  9k,j±1  =  0  • 


(See  Figure  2).  This  result  follows  immediately  from  the  fact 


corner  points  of  lie  in  RD. 

D 


s  useful  to  write  L^^  as  the  sum  of  two  operators,  one  essentially 
and  the  other  on 


H  2h 


4 
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Definition: 


(3.24a) 

(3.24b) 

(3.25a) 

(3.25b) 


Lemma  3.5; 


(3.26) 


Let  M  ,  Q  :  Su  -*■  Su  be  defined  by 

X  X  n  n 


tQxU^k,j  =  0  ’  (xk,yj^  £  Vn2h  ’ 

^xU^k,j  =  ^'LH  ^k,j  5  (xk,yj)  e  fi2h  * 

tMxU^k,j  =  *  (xk’yjJ  £  Vfi2h  ’ 


*Wj  =  0  »  (xk,yjJ  £  Q2h  ‘ 


Let  v  e  S2h  .  Then 


i<V2h,-I2h,>H  *  <L2hv'v>2hl  lKh2<L2hv-v> 


Proof:  The  lemma  follows  from  Lemma  3.4  and  the  estimates 

|4*k+l,j  '  pk+l,jl  -  Kh2pk+l,j  ’ 

1 4®k, j+1  ‘  pk,j+l 1  -  Kh2pk,j+1  • 


Remark:  When  p(x,y)  =  1  ,  then  K  =  0  . 

Finally,  we  "lift"  L2h  (an  operator  defined  on  S2h)  as  follows: 

L  :  S0.  +  Su  be  defined  by 
x  2h  H  J 

(3.27a)  ^x^2hv^kj  =  °»  ^xk,yj^  6  QH/f22h  ’ 

(3.27b)  ^x(l2hv^k,j  =  ^L2hv^k,j  *  (xh’yj)  €  Q2h  ‘ 

Remark:  Using  this  definition  we  may  rephrase  (3.26)  as 

(3.28)  I {  QxI2hV,I2hV>H  "  ( '■xI2hv,I2hV>H^  -  Kh  (  ‘‘xI2hV  ,!2hv  ^  • 


1  et 


4.  Analysis  I. 


Let  =  u  -  be  the  initial  error.  Then  e  =  u  -  u  is  the  error 
after  step  1,  the  smoothing  step.  Assumption  (2.16)  asserts  that 

(4.1)  II  e  11  l  =  (  Lhe’£>  1  <  L^.e0)  =  ||  e°||^  . 

h  h 

Using  the  decomposition  (3.4a)  we  have 

(4.2)  e  =  nh  +  w  ,  nh  e  Nh  ,  w  €  S|_|  . 

From  step  2  [i.e.  (2.18)]  of  the  algorithm  and  Lemma  3.1  [i.e.  (3.5b)]  we 
see  that 

(4.3)  e  =  u-  u  =  l|jw. 

Hence,  using  (3.4a)  we  see  that 


Using  (4.3)  and  (3.7a)  and  step  3  of  the  algorithm  we  see  that 

(4-5)  v  ■  (IM>“  ■  rH  • 

See  [3]  for  a  more  complete  discussion  of  the  significance  of  this  fact. 
Lemma  4.1:  Let  v  e  SH  be  the  solution  of 

A 

(4.6)  Lhv  =  rH  =  Lhw  . 

Let 


(4.7) 


v  *  "H  +  I2hV  ' 


nu  e  Nu  , 


V  e  S 


Let  i p  be  the  function  in  SH  constructed  in  step  4  [i.e.  (2.19)]  of 
the  algorithm.  Then 


(4.8) 


ip  = 


nH  • 


Proof:  Observe  that  (2.19a)  and  (3.5b)  imply  that  $  e  N, ,  .  Also  (2.19b) 
and  (4.6)  yield 

[LH(v-ij))]k^  =  0  ,  (xk>Yj)  e  . 

That  is 

(4.9a)  (v-40  =  [ (nH“4>)  +  1 2hV-^  £  KH 

whi  1  e 

(4.9b)  (nH-<£)  e  Nh  . 

Using  (3.4a)  and  (3.4b)  we  see  that  (4.8)  holds.  I 


Consider  the  function  <J>  which  is  constructed  in  step  6  of  the  algorithm 
We  have 

(4.10)  l2h0  =  IH  LH(v-$)  =  IH  LHI2hV  • 

thus 


(4.11) 


"2fr 


V  . 


From  (4.3),  (4.11)  and  step  7  of  the  algorithm  we  see  that 
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(4.12) 


e1  =  u  -  u1  =  I^[(w-$)  -I^]  e  Kh  * 


0 


Thus,  if  we  seek  an  eigenfunction  e  ,  it  must  have  the  form 


0  Th 
e  -  IHeH  . 


As  we  shall  see,  the  generality  of  and  the  estimate  (4, 
that  it  suffices  to  consider  the  case  where  =  1^  .  In  that 


(4.13)  l  =  e°  =  lj][nH  +  1^11]  ;  nu  e  nh  ,  u  £  s?h  . 


'H  H 


2h 


If  e1  =  ye0  (4.12)  becomes 


e1  =  i|][(nH-$)  +l2h(u-<J>)3  =  lil|jCnH  +I2hU^  * 


Thus 


(4.14) 


4>  =  AnH  ,  <P  =  Au  ,  A  *  (1-p) 


Returning  to  Lemma  4.1  we  have 


(4.15) 


lh<*+4v>  ■ 


LH(AnH  + i2hv) 


hl«H  +I2hU) 


LH^H  +I2hU) 


From  (3.8b),  (3.6a),  (4.11)  and  (4.14)  we  see  that 


(4.16a) 

(4.16b) 


LHI2hV|0  -  24hV  *  2Lh*  ’  2XL2hU 
u2h 


L  v 
LH12hv 


=  0  . 


VQ2h 


1)  implies 
case 


Thus,  (4.16)  and  the  definition  of  L^Li.e.  ( 3 . 27 ) J  allows  us  to  rewrite 
(4.15)  as 

(4.17)  X[LHriH +2LxI2hU]  =  LH^H+I2hU^  • 

To  simplify  the  eigenvalue  problem  (4.17)  we  define 

.  #  c  c 
L  :  SH  ^  SH 

as  follows:  let  v  e  Su  .  Then  there  is  a  unique  representation 

H 

(4.18a)  v  =  CH  +  I^W  ,  CH  e  Nh  ,  W  «  S2h  . 

Then 

(4.18b)  L#v  =  Lh?h  +  2LxI2hW  . 

The  eigenvalue  problem  (4.18)  now  becomes 

(4.19a)  AL#v  =  L^v  , 

(4.19b)  v  =  nH  +  J2hU  • 

#  A 

Observe  that  both  L  and  LH  are  symmetric  positive  definite  operators 
Therefore,  there  is  a  complete  set  of  eigenfunctions  {v^}  which  satisfy 

(4.20)  <L#vk,Vj>  =  <  L^v^,Vj  >  =  0  ,  M  J  • 

Then  (3.7c)  implies  that 


(4.21) 


£  max  | 1 -A |  =  max  |y|  . 
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Thus,  in  view  of  (4.1),  the  general  three-grid  iteration  (Gh  f  Ih) 
also  satisfies  (4.21 ) . 
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5.  Analysis  II 

Consider  the  basic  eigenvalue  problem  (4.19).  Let  us  now  focus  our 
attention  on  the  left-hand-side  of  (4.19a).  Using  (3.9a)  and  (4.19b)  we 
have 

(5.1a)  .  1  L<'>Sh  *  1  l'1^  *  1  |  hS’)l2hU  • 

and 

<v,Lhv>  =  2  (  VLH  V  +  2  ( VLH  ^2hU>  +  2  *  V*>i  > 

(S.lb)  +  2  <^H’LH  ^2hU>  +  2  { I2hU,LH  \  }  +  2  (  I2hU,LH  ^2hU  > 

1  /  .H  „  r(l)-  N  A  1  /  tH  ..  r (1  )tH  ,.v 

2  I2hU*LH  nH  2  I2hU,LH  !2hU  ’ 

The  basic  estimate  (3.12a)  allows  us  to  replace  by  LH  provided 

we  accept  error  terms  of  the  form 

(5.2a)  =  Kh2[<LHl2hU,I2hU>  <LHnH,nH>]l5> 

(5.2b)  fi2  =  Kh2<  LHI2hU,I2hU>  ’ 

(5.2c)  63  =  Kh2  ( lHnH,nH>  • 

Thus  we  may  rewrite  (5.1b)  as 

<'-V)4<SH'LHfiH>*7<I2hU'LHI2hU> 


*  2  ^H’^H  >  +  ^H’^H  ^ I2h°  >  +  2  (  I2hU,‘’H  *!2tiU>  ' 


0(6) 


(5.3) 


5.4) 


0(6)  =  0(61  +52  +63)  . 
From  (3.6b)  of  Lemma  (3.1)  we  see  that 


(Vk,j  =  0  * 

<vV  6  1 

ience 

(5.5a) 

(V^H  ^H> 

<S„,mxSh>  . 

(5.5b) 

<V"»IaM> 

Thus,  we  may  rewrite  (5.3)  as 


<  v.CHv  > .  1  <  vLhSh  >  ♦  1  <  iV.LHi"hu>  +  \  <S„.MXSH  > 

(5.6)  ♦<VV»U>  +  I<4U-Va.U>  + 

T<I2hU-V2l,U>  +  0({)- 

Let  us  consider  the  term 

(5.7a)  J:  *  \  ( I2hU,LHI2hU ' 

From  (3.7b),  (2.12)  and  (3.9b),  (3.19b)  we  have 

(5.7b)  J  =  <  U,L2hU  >2h  -  j  <  U’L2hU  )  2h+  2  (  U’^2h^U>2h  ‘ 
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Thus,  using  the  definition  of  Lx  and  (3.17c)  we  obtain 


(5.7c) 


J  =  2  (  I2hU,^xI2hU  }H  +  2  <U’^2h^U>2h 


The  estimate  (3.28)  allows  us  to  replace  Q  by  L  provided  we  accept 

X  X 

errors  of  the  form 


(5.8)  5  =  Kh2  < LxI2hU’I2hU >  =  Kh  (  LHI2hU,I2hU>  ’ 


Thus,  we  rewrite  (5.6)  as 


(5.9) 


<  v,Lhv  >  -  j  <  v,L#v  >  +  j  <  U.L^U  >2h 


+  j  <v,Mxv>  +  0(6)  +  0(6)  . 


The  eigenvalue  problem  (4.19)  becomes 


(5.10) 


U-Jj)  <  v,L#v>  =  \  <U,L^)U>2h  +  \  <v,Mxv>  +  0(6+6) 


Hence 


X-l>  -Kh2 


(5.11) 


A  >  j  (1-Kh2) 


The  complete  proof  of  our  basic  estimate  requires  a  more  detailed 
analysis  of  the  terms  which  appear  in  (5.10). 


V  .*•  .*»  „*v  A  *,  *. 
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Lemma  5.1:  For  all  n  e  NH  we  have 

(5.12)  0  <  < Mxn,n > £  (1+Kh)  < LHn»n>  . 

Proof:  Using  the  description  of  given  in  Lemma  3.3  we  see  that 

(5.13)  <Mxn,n><  l  pkj.(l+Kh)(nkj.)2  . 

On  the  other  hand,  since  nkj-  =  0  if  (x^.y^)  e  ,  (2.13b)  shows  that 

(5.14)  <LHn,n>  =  ^2  I  ckj(nkj)2  -  £  pkj (1 'Kh2) (nkj )2  * 

Thus,  the  lemma  is  proven. 

Our  next  result  is  intuitively  clear.  Nevertheless,  the  details  are 
somewhat  technical. 

Lemma  5.2:  For  every  U  e  ,  we  have 

(5.15)  <4U'Nx,H2hU>i'UK,’l(I"hU'Cx4U>  ■ 


Proof :  The  idea  behind  the  proof  is  quite  simple.  We  have 


fLxI2hU^kj  =  £4hUVj 


(xk,yj)  e  n2h 


while  [see  (2.11)]  {(I^U^  for  ^xk’yj^  6  ^H^2h  is  an  "avera9e"  of 

(U  }  with  (x  ,y  )  e  .  Since  M  and  L  are  "almost"  the  same 


operator,  (5.15)  should  follow.  The  complete  details  of  the  proof  are 
given  in  the  Appendix.  | 

Theorem  5.1 :  Consider  the  three  grid  iterative  scheme  described  in 
section  2:  steps  1-8.  Let 

0  0  1  1 
£=U-U,  £  =  U  -  U 

There  is  a  constant  K  _>  0  ,  depending  only  on  p(x,y)  and  its  first 
and  second  derivatives,  such  that 

(5.16)  ||  e1  ltL  0+Kh)||  e°||L  . 

h  h 

Moreover,  if  p(x,y)  =  const  >  0  ,  then  K  =  0  . 

Proof:  Let  (A,v)  be  an  eigenvalue  and  eigenfunction  of  (4.19),  or 
equivalently,  (4.17).  As  we  have  seen,  (A,v)  satisfy  (5.10)  and  (5.11) 
holds.  Expanding  the  terms  of  <v,Mxv>  we  have 

(5.17)  R-H.S  =l<U,l'2)u>2h  +i[<SH.MxSH>  +2(SH,Mxl"hU>  *<IH2hU,Hxl”hU>] 

+  0(6+6)  . 

Using  (3.19b),  (3.27),  Lemma  5.1  and  Lemma  5.2  we  have 

R.H.S  <  (1+Kh)  <  I^U,LxI^U)  +  <nH,MxnH>  +  <  l2hU’MxI2hU>  +  0(6+^ 

<  (1+Kh)  [2  <  I2hU’^xI2hU  J  +  (  VLhV]  +  0(6+^  * 
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Thus 


(A  -  %)  <  v,L#v  >  <  0(6+6 )  +  Kh  <  v,L#v>  . 


Hence 


(5.18)  A  <  |  (1+Kh)  . 

Thus,  (5.18)  and  (5.11)  together  with  the  remarks  at  the.  end  of  section  4 
imply  the  theorem. 
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Appendix 


In  this  appendix  we  give  the  details  of  the  proof  of  Lemma  5.2.  We 
use  the  formulae  of  Lemma  3.3  which  give  the  form  of  L^  and  hence 
give  the  form  of  .  A  simple  summation-by-parts  argument  shows  that 


( A  *  1 )  <  I2hU,MxI2hU>  =  h2  2  *k+l  ,j^I2hU)k+2,j  "  (I2hU)k,j] 

G 

+  72  ^  Vj+l£(I2hU)k,j+2  "  (I2hU)k,j-*  ’ 

n  kg 


whi  1  e 


(A. 2) 


<  1^  U  L  1^  U  > 
1 2hU ’  x  2hU ' 


+ 


n  Pk+1  ,  j 

RB 

n  Pk,j+1 

kb 


[U 


k+2,j 


[U 


k,  j+2 


■Uk,j 


-U.  . 
k,j 


Note:  It  is  essential  to  observe  that  the  sums  in  (A.l)  are  taken  over 


points  (x. ,y.)  e  R.  while  the  sums  in  (A. 2)  are  taken  over  points  in 
K  J  b 


*B  .  Moreover,  if  (x^y^)  e  Rg 


then 


(,2hU)kJ  *  Vj  • 


Let  j  =  1  (mod  2)  be  fixed  and  consider  the  contribution  to  the  first 
sum  (on  the  right-hand-side)  of  (A.l)  from  the  points  on  a  connected  segment 
of  the  intersection  of  ft  with  the  line  y  =  jh  .  That  is,  we  consider  a  sum 
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(A. 3) 


h  =4  i 


r=r. 


A2r,jt(I2hU)2r+l  ,j 


(I2hU)2r-l 


For 


rQ  £  r  <  r,  ■  1  all  the  points  (x2r+l’y'^  €  ^  *  Ttie  P01nts 


(x 


2r  .T»yj)»  (x2r  +-j  »y  j )  may  lie  on  3ft  or  may  lie  outside  ft 


n(*) 


depending  on  the  slope  of  the  boundary  near  these  points.  In  either 

CaSe  (I2hU>2Vl,j  '  «»%♦!, J  -°- 

Consider  the  contribution  of  such  an  "end-point"  to  the  sum 
(See  Fig.  3).  For  definiteness,  consider  the  term 

(A>4)  J2rQ,j  =  A2r0,jC(I2hU)2r0+l,j  “ (I2hU)2rQ-l , j]2  * 

We  should  consider  two  cases,  either  the  boundary  has  slope  °°  near 

(x„  y.)  or  slope  +  1  .  When  the  slope  is  +  1  ,  then  6,  .  =  0  and 

zr0’  J  ~  "  Zr0’J 

the  analysis  is  much  like  the  case  when  both  (x2r+i’yj)»  (x2r  l’yj^  e  n  • 

However,  when  the  slope  is  °° ,  then  0_  .  ^  0  (see  Fig.  2).  In  this 

^rQ,j 

case 

(I2hU>2rn-T,d  ">• 


l(I2hU)2rntl,jl  -  ^  4 hl  C|U2r„+2,j+l'  *  lU2rn+2,j-l I ] 


and 


0  -  A2rQ,j  -  8  (1+Kh)p2r0+l,jil 


Near  a  reentrant  corner  with  angle  45°  one  of  these  points  may  actually 

lie  on  3ft  while  the  segment  from  its  neighbor,  say  (x7  ,  ,y.),  is 

_  zr0  J 
not  entirely  in  ft  .  Nevertheless,  U  =  0  at  that  point. 


Totally  interior 


Near  the  boundary 


Near  an  oblique  boundary  edge 

■  denotes  a  point  in  Rg 
□  denotes  a  point  in  R 


Figure  3 
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Hence 


J2rn,j  —  16*4  (1+Kh)(P2rn+l,j+l[U2rn+2,j+l  " U2rn. j+1  -1 


(A. 5) 


+  p2rn+2,j-l[U2rn+2,j-l  'U2rn,j-1] 


Now,  consider  a  term 


(A-6)  J2r, j  =  A2r,jt(I2hU)2r+l,j '(I2hU)2r-1,j-1 


where  rg  <  r  <  r]  •  We  wrlte 


!r,j  =  (I2hU)2r+l,j  "  (I2hU)2r-l,j 


as  a  sum  of  four  terms 


I  .  =  D,  ,  ,  +  D0  .  ,  +  Di  .  ,  +  D0  .  , 
r,j  1 ,J+1  2,j+l  1  ,  j-1  2,j-l 


where  -  up  to  terms  of  order  h‘ 


fa  7,  \  n  -  ^2r+3/2 ,  J+h  1 1  _  <•>  •-?  y 

(A'7a)  4p2r+1J  'Wj.l  4P2r.,_.  U2r,j+1 


P2r-^j_^  j. 


n  -  ^r+Sj.j+h  P2r-3/2,j+j|  y 

2'J+1  4P2r+l ,  j  2r'J+1  4p2r-l,j  2r'2’J+1 


(A. 7b) 
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Indeed,  expanding  the  coefficients  about  (x^.y.)  yields 


V.jl  -  1  4^  {'U2r+2,j+l  -U2r,j+ll  +  '  U2r ,  j+1  'U2r-2,j+l 


+  lu,  ^  .  ,  -U,  .  ,|  +  |U0  .  ,  -U,  0  .  , |}  . 

1  2r+2, j-1  2r,j-l  1  2r,j-l  2r-2,j-l 1 


Hence 


J2r,j  -  MT  {p2r+l,j+l[U2r+2,j+l  _U2r,j+l]2  + 


p2r-1,j+l[U2r,j+l  'U2r-2,j+l]2  +  P2r+1 , j-1 [U2r+2 , j-1  "U2r,j-1]2  + 


P2r-1,  j-1 [U2r, j-1  ' U2r-2, j-1 ]2}  ' 


Upon  adding  the  contribution  from  each  j  line,  we  see  that  each 

term 

1  2 
^2  P2r+l,j+l[U2r+2,j+l  'U2r,j+1] 

enters  at  most  4  times.  Since  each  such  term  has  a  coefficient  which  is 
less  than  or  equal  to  (l+Kh)/4,  the  lemma  is  proven.  | 
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